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^ ■ 1 Introduction 

\ Since the 1964 discovery of K^ — > -k + tt~ [[!]], understanding CP violation has been a 

major goal of particle physics. Subsequent experiments showed that the dominant 
mechanism of CP violation in neutral kaons (still the only particle system in which CP 



violation has been observed) is a small asymmetry between K° — > K and ~K° — > K° 
mixing. This asymmetry, referred to as indirect CP violation, results in the Kl and 
^j- . Kg being states of mixed CP. The parameter e, which is used to parametrize this 

effect, quantifies the CP impurity of the Kl and K$ states: 

Kl ~ -ft'odd + tK even , 

X 

<D ! where |e| = 2.28 x 10 ^ CP\K even >= +l\K even >, and CP\K odd >= -l|if odd > . 

o 1 There has been considerable effort during the last 30 years to determine whether 

or not the CP symmetry is also violated in the decay amplitude (i.e., K odd — > nn). 
This effect is referred to as direct CP violation and is parametrized by e'. The ratio 
e'/e can be determined from the double ratio of the 2-pion decay rates of Kl and Ks'. 

% : Re(e7e) „ I \ nx L -» *+*-)/mc 8 -» tt^q _ 



(1) 



e'/e 7^ is an unambiguous indication of direct CP violation. 

The Standard Model predicts both direct and indirect CP violation. Unfortu- 
nately, large hadronic uncertainties make precise calculations of Re(e'/e) difficult. 
Most recent Standard Model predictions |2j are in the range Re(e'/e) = (0 — 30) x 10 -4 . 
Other models, such as the Superweak Model of Wolfenstein ||, predict no direct CP 
violating effects. 

The two best previous measurements of e'/e come from E731 0] at Fermilab and 
NA31 [g at CERN: 

Re(e'/e) = (7.4 ± 5.9) x 10~ 4 (E731) 
Re(e'/e) = (23 ± 6.5) x 10" 4 (NA31). 
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The CERN result is 3.5 standard deviations from zero, while the Fermilab result is 
just over 1 sigma from zero. To clarify this experimental situation and definitively 
resolve the question of whether or not direct CP violation occurs, new experiments 
were built at Fermilab (KTeV), CERN (NA48), and Frascati (KLOE) to attempt 
to measure e'/e at the (1 — 2) x 10 -4 level. The KTeV and NA48 experiments are 
similar fixed-target experiments. They differ mainly in the method used to produce 
K$, and in the technique used to correct for the difference in detector acceptance for 
K$ and decays resulting from the large Ka, — lifetime difference. The KLOE 
experiment at Frascati is trying a new technique using an e + e~ — > collider. 

The first measurement of e'/e from the KTeV experiment || is the subject of this 
paper. A more complete description of the KTeV measurement is given in |7]]. The 
NA48 and KLOE experiments are described by Barr || and Bertolucci |J in these 
proceedings. 

2 The KTeV Detector 

To achieve the required level of statistical and systematic uncertainty in e'/e, the 
KTeV experiment (Fig. |l|) uses the same double-beam technique as E731 with a 
new detector and beamline. Following the primary target, collimators and sweeping 
magnets are used to form two almost parallel neutral beams. A fully active regenerator 
is placed in one of the beams 122m from the production target, at the upstream 
end of the decay region, to provide a source of K$ for the experiment (this beam 
is referred to as the regenerator beam and the other as the vacuum beam). The 
regenerator (along with a movable absorber that attenuates the beam hitting the 
regenerator) is moved from one beam to the other each minute to eliminate many 
possible systematic errors in normalization and detector response. All four K — > 2n 
decays are detected simultaneously. The detector consists of a large vacuum decay 
region instrumented with photon veto counters, a drift chamber spectrometer, and 
a Csl electromagnetic calorimeter. Compared to E731, KTeV also has an improved 
trigger and data acquisition system. The final stage of the trigger includes full event 
reconstruction and filtering before data are written to tape. 

The performance of the calorimeter, made up of 3100 pure Csl crystals, is particu- 
larly important to the success of the experiment. The energy scale of the calorimeter 
is directly related to the reconstructed decay position along the beam (z) direction 
for K — > 7T 7r° decays, and is therefore a critical part of understanding the detector's 
acceptance for neutral events. Figure shows the energy resolution of the calorimeter 
as a function of momentum for electrons from K — > irev events. The average energy 
resolution for photons from K — > tt tt° events is 0.7%. The excellent energy resolution 
also reduces background for both the 7i + ii~ and 7r°7r° decay modes. 
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Figure 1: Diagram of the KTeV detector. 



3 K — > 7t + tt and K — » n tt Analysis 

The analysis presented here is based on K — > 7r°7r° data collected during 1996 and 
K — > 7r + 7r~ data collected during 1997 fllQf . Using charged and neutral data from 



different running periods does not significantly increase the systematic error in e'/e 
because the two decay modes use essentially independent detector systems. It is 
critical, however, that the Kg/K^ flux ratio be the same in both years. This issue 
will be addressed in Section 4. 

K — > 7r + 7r~ events are reconstructed in the charged spectrometer, consisting of 
four drift chambers, two on either side of a large dipole magnet providing a 0.41 
GeV/c horizontal momentum kick. Important requirements for the 7r + 7r _ decay mode 
include: 

• each track must have a momentum of at least 8 GeV/c and deposit less than 
85% of its energy in the Csl calorimeter; 

• there must not be any hits in the muon hodoscope located behind 4 m of steel; 

• the square of the transverse momentum of the 7r + 7r~ system relative to the 
initial kaon trajectory, p 2 , is required to be less than 250 MeV 2 /c 2 . 
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Figure 2: Csl energy resolution as a function of momentum measured with 200 
million electrons from K —>■ Treu decays. 

The K -> vrV reconstruction is based on the energies and positions of four 
photons measured in the Csl calorimeter. The events are reconstructed by selecting 
the photon pairing which is most consistent with both 7r° decays occurring at the same 
point. As an alternative to p 2 , which cannot be reconstructed in 7r°7r° decays because 
the photon directions are not measured, we calculate a "ring number" based on the 
center-of-energy of the four photons at the calorimeter. The ring number is defined 
as the area (in cm 2 ) of the smallest square that is centered on the beam and contains 
the center-of-energy. We require that the ring number be less than 110, which selects 
events with center-of-energy inside a square region of area 110 cm 2 centered on each 
beam. 

Invariant mass distributions for the K — > 2n decay modes for events with 110 < 
z < 158 m and 40 < Er < 160 GeV are shown in Fig. |^. The corresponding 
distributions of decay positions along the beam (z) direction are shown in Fig. 

There are two classes of background in these K — > tttt samples: misidentified 
kaon decays and real K — > -kit events that have scattered in the regenerator or final 
collimator. The total background levels for the four decay modes are shown in Fig. || 
Table | summarizes the different components of the background. For the tt + tt decay 
mode, backgrounds in both beams are less than 0.1%. In the vacuum beam, the 
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Figure 3: K — > txtx invariant- mass plots before background subtraction. 
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Figure 4: Decay vertex distributions for (a) K — > 7t + tt and (b) K — > 7r 7r decays, 
showing the difference between the "regenerator" (K$) and "vacuum" (K^) beams. 



background comes mainly from misidentified semileptonic decays. In the regenerator 
beam, the main background is from kaons that scatter in the regenerator before 
decaying to 7r + 7r~. The background levels are much higher for the 2ir° decay mode. 
As in the charged decay mode, kaons that scatter in the regenerator are the main 
background in the regenerator beam. Since the ring-number variable is not as effective 
as pi at identifying scattered kaons, however, the neutral-mode background from this 
source is 1.07% (more than an order of magnitude larger than in the charged decay 
mode). Kaons that scatter enough to be reconstructed in the wrong beam contribute 
a background of 0.3% to neutral decays in the vacuum beam. The vacuum beam 
also has a background of 0.27% from — > 3ir° decays with lost and/or overlapping 
photons. The numbers of K — > 2tt events after background subtraction are given in 
Table |. 
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Table 1: Background levels (in %) for K — > 7i + n and K — > tt tt' 
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Figure 5: Distributions of pf for the 7r + 7r~ samples and ring number for the 7r°7r° 
samples. Total background levels and uncertainties are given for the samples passing 
the analysis requirements indicated with arrows. 
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Vacuum Beam (Ki,) 


Regenerator Beam (Kg) 


K — > 7T + 7T 
K -> 7T°7V° 


2,607,274 
862,254 


4,515,928 
1,433,923 



Table 2: Numbers of events after background subtraction 



As shown in Fig. |4], the difference between the and K$ lifetimes results in 
very different decay vertex distributions for the and K$ decays which must be 
compared to compute e'/e. Therefore, to extract a value for e'/e, the numbers in 
Table |2] must be corrected for the variation in detector acceptance as a function of 
z. We make this correction with a Monte Carlo (MC) simulation. The simulation 
models kaon production and regeneration to generate decays with the same energy 
and z distribution as the data. It also includes a detailed simulation of all detector 
elements. 

The quality of the Monte Carlo simulation is studied using distributions from 
K — > 2tt decays, as well as higher statistics decay modes. Figure |^ shows a comparison 
of data and Monte Carlo vacuum beam z vertex distributions for the tttc signal modes, 
as well as for the much larger nev and 3tt° samples. Since the average decay positions 
for i^L and K$ decays differ by about 6 m, a relative slope of 10~ 4 per meter in the 
data/MC ratio would result in an error of 10~ 4 in e'/ e. The only noticeable problem in 
Fig. | is the slope of (-1.60 ±0.63) x 10" 4 m" 1 in the data/MC slope for K -> tt+tt-. 
Although this slope is only 2.5 sigma from zero and the slope in Kl — > nev is much 
smaller, we assign a systematic error of 1.6 x 10~ 4 on Re (e'/e) based on the full size 
of the observed slope in Kj, — > 7r + 7r _ . The 7r°7r° and 37r° data and Monte Carlo 
z distributions are consistent. Since the 37r° decay mode is more sensitive to most 
acceptance problems, we use the slope in the data/MC ratio in 3ti to place a limit 
of 0.7 ± 10~ 4 on the bias in Re(e'/e) from the neutral-mode acceptance. 

4 Systematic Errors 

Table |3| summarizes the estimated systematic uncertainties in e'/e. The 7r + 7r~ and 
7r°7r° columns list errors that are specific to the charged and neutral decay modes. 
The errors at the bottom of the table are common to both decay modes. Adding 
these errors in quadrature yields a total systematic error on e'/e of 2.8 x 10~ 4 . Only 
a few of these errors will be addressed here; additional details are given in @. 

The largest contribution to the systematic error comes from uncertainties in the 
z dependence of the acceptance, which are estimated from the comparison of z dis- 
tributions for data and Monte Carlo discussed above. Other significant uncertainties 
result from the energy scale and background subtraction in neutral mode. To eval- 
uate the uncertainty in the neutral energy scale, we compare the reconstructed z 
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Figure 6: (a) Data - Monte Carlo comparisons of vacuum-beam z distributions for 
7r + 7r~, neu, 7r°7r°, and 37r° decays; (b) straight-line fits to the data/MC ratios. 
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Table 3: Summary of systematic 


uncertainties 


in 


Re(e'/e) 



vertex of 7r°7r° events produced at the downstream edge of the regenerator with the 
reconstructed z vertex for 7r° pairs produced by hadronic interactions in the vacuum 
window at the downstream end of the decay region. The difference between these 
measurements is 2 cm greater than the actual distance between the regenerator edge 
and the vacuum window, resulting in a systematic error of 0.7 x 10~ 4 on Re(e'/e). 
The neutral mode background uncertainty results largely from uncertainty in the ac- 
ceptance for scattered K — > 7r + 7r~ decays, which are used to tune the Monte Carlo 
simulation for scattered events. 

As mentioned earlier, in combining 7r + 7r~ and 7r°7r° data from different years, we 
must consider the possibility of a change in the Kg/K^ flux ratio between the two 
samples. Although the same regenerator and movable absorber were used for the two 
data samples, we assign a small uncertainty on Re(e'/e) to account for the possible 
effect of a temperature difference between the two data collection periods, which 
would change the densities of the movable absorber and regenerator. 

5 Results 

The numbers of events and relative acceptances (from the Monte Carlo simulation) 
for the four 2tt decay modes can be used to calculate a simple estimate of e'/e using 
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Figure 7: Decay vertex distribution from K — > 7t + tt~ events with 30 < Ek < 35 GeV 
downstream of the regenerator showing K$ — interference. 



Equation 1. This calculation JL1]] yields Re(e'/e) = (26.5 ± 3.0) x 10 , but is not 



precisely correct because it assumes that there are only Kg decays in the regenerator 
beam. As illustrated in Fig. [7], the regenerator beam contains a coherent mixture of 
Kg and K^, which must be taken into account in the calculation of e'/e. The decay 
distibution in the regenerator beam also allows us to measure the K$ lifetime (js), 
the K$ — mass difference (Am) , and the relative phases of the CP violating and 
CP conserving amplitudes ($ o f° r K — > 7r°7r° and $ ++ for K — > 7r + 7r~): 

t s = (0.8967 ± 0.0007) x 10~ 10 s 
Am = (0.5280 ± 0.0013) x lO^^s" 1 
A$ = $oo - = 0.09° ± 0.46°. 

The final value of Re(e'/e) is extracted from the background-subtracted and 
acceptance-corrected data with a fitting program that calculates decay vertex dis- 
tributions, properly treating regeneration and K$ — interference downstream of 
the regenerator. Including the systematic error from Table [3], the result of the fit is 

Re(e'/e) = (28.0 ± 3.0 (stat) ± 2.8 (syst)) x 10~ 4 
= (28.0 ± 4.1) x 10~ 4 . 



We have performed several cross checks on this result. Consistent results are 
obtained at all kaon energies, beam intensities, periods during the run, magnet polar- 
ities, and for both regenerator positions. We also have done the analysis using 7r + 7r~ 
data from 1996 (collected simultaneously with the 7r°7r° data) rather than 7r + 7r~ data 
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Figure 8: Re(e'/e) as a function of kaon energy. The horizontal line represents the 
central value of Re(e'/e) = 28.0 x 10 -4 . Only statistical errors are shown. 



from 1997. The result is consistent with that quoted above, but with a larger sys- 
tematic error [10]. Some of these cross checks are summarized in Figs. |8| and |9|. It 
is worth noting that the e'/e analysis was done "blind": the value of e'/e was hidden 
with an unknown offset until after the analysis and evaluation of the systematic error 
were completed. 



6 Conclusions 

Based on about 1/4 of the data collected during the Fermilab 1996-1997 fixed-target 
run, KTeV has measured Re(e'/e) = (28.0±3.0 (stat)±2.8 (syst)) x 10" 4 . This result 
establishes the existence of CP violation in a decay process at almost 7 sigma, and 
rules out the Superweak Model as the sole source of CP violation. Although the result 
is larger than most Standard Model calculations, it supports the Standard Model 
explanation of CP violation. Figure [10] shows a comparison of the KTeV measurement 
with previous results; the preliminary result of the NA48 experiment || O] is shown 
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also. The KTeV result is consistent with earlier evidence for direct CP violation from 
NA31 and differs from the E731 result by 2.9 sigma. Study of the E731 result has 
not revealed any error that would cause this discrepancy. A weighted average of all 
e'/e measurements gives Re(e'/e) = (21.3 ±2.8) x 10~ 4 with a confidence level of 7%. 

The analysis of the remaining 3/4 of KTeV's 1996-1997 data sample is in progress. 
In 1999, we collected a new e'/e data set with statistics equal to the full 1996-1997 
sample. Several small detector modifications were made to improve the systematic 
quality of the 1999 data. We also performed additional systematic studies during the 
1999 run. The full data sample should allow us to reduce the statistical uncertainty 
on Re(e'/e) to 1 x 10 -4 . Significant work will be required to reduce the systematic 
error to a similar level. 

Although the very existence of direct CP violation has been an issue until recently, 
the full data sets of KTeV, NA48, and KLOE may soon provide measurements of e'/e 
at the 5% level. Considerable improvement in theoretical calculations of e'/e will be 
required to take full advantage of this precision. There is some optimism, however, 
that the next rounds of calculations using lattice gauge theory may approach a 10% 
uncertainty in e'/e, making the precise measurements of e'/e equally precise tests of 
the Standard Model. 
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Discussion 



Sherwood Parker (University of Hawaii): Did you make any changes to the 
apparatus to achieve the better results in the second run? 

Blucher: There were no fundamental changes to the detector, but small modifica- 
tions were made to improve data quality and datataking efficiency For data quality 
the most important changes involved the drift chamber system. The drift chamber 
electronics were improved to have higher gain and lower noise to reduce the chamber 
inefficiency. The two upstream drift chambers were also restrung because of damage 
sustained during the 1996-1997 run. 

B. F. L. Ward (University of Tennessee): Is there any understanding of why 
there may have been a systematic bias in the earlier measurement that caused it to 
miss the higher value of Re(e'/e)? 

Blucher: As mentioned earler, we have studied the E731 analysis and have not 
found any evidence of a systematic problem. For example, comparisons of data and 
Monte Carlo z distributions for E731 [T3|], similar to those shown for KTeV in Fig. 
do not show any sign of an acceptance problem. 

Mario Calvetti (INFN, Florence): In your analysis, you rely on the Monte 
Carlo acceptance corrections to the first order. Why is this preferable to the method 
of NA48; that is, reweighting the events in order to have similar longitudinal vertex 
distribution? 

Blucher: The NA48 procedure sacrifices statistics to reduce the required accep- 
tance correction, while the KTeV procedure maximizes the statistical power of the 
data sample, but requires that the detector acceptance be understood. We believe 
that we have a reliable procedure to estimate the systematic error associated with 
our acceptance correction. As a cross check, we have analyzed our data with an alter- 
nate technique that compares vacuum and regenerator beam z distributions directly, 
eliminating the need for a Monte Carlo acceptance correction. This analysis gives a 
consistent value of e'/ e but with a significantly larger statistical error. 
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